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We interpret the KARMEN time anomaly as being due to the production of a (dominantly bino) 
neutralino with mass 33.9 MeV, which is the lightest supersymmetric particle but decays into 3 
leptons through the violation of R-parity. For independent gaugino masses Mi and Ms we find 
regions in the (Mi, M2, fi, tan/3) parameter space where such a light neutralino is consistent with 
all experiments. Future tests of this hypothesis are outlined. 

o 
o 
o 

I. INTRODUCTION 

o . 

["T 1 , In 1995, the KARMEN experiment at the Rutherford Appleton Laboratory reported an anomaly in the time 
distribution of the charged and neutral current events induced by neutrinos from 7r + and n + decays at rest [Q . This 
was ascribed to the production of a new particle, denoted 2, in the anomalous pion decay 

7T+ (1) 

• with a small branching ratio in the range ~ 10 — 10 depending on the lifetime of x. The particle must be neutral 
C<~) ' since it passes through over 7 m of steel shielding. The time-of-fiight to the detector is 3.6 ±0.25 fis, implying that the 
particle moves non-relativistically with velocity v x = 5. x 10 6 ms" 1 . This requires its mass to be 33.9 MeV and 
its kinetic energy to be T x w 5 keVQ The observed energy in the detector is ~ 11 — 35 MeV, which must therefore 
come from the decay of x. Since 1995 the KARMEN experiment has been upgraded to significantly reduce the cosmic 
ray background. It has recently been reported that the time anomaly persists in the new data A time-of-flight 
likelihood analysis adopting the hypothesis that it is due to a decaying particle as described above has a negative 
natural log-likelihood ratio of 9, i.e. less than 1 in 10 4 chance of being a statistical fluctuation. The significance is 
f"*> || thus sufficiently high that we are motivated to reexamine its physical origin. 

There have already been several proposals to explain the KARMEN anomaly (| ||. In Rcf. j|, the authors 
considered in detail the possibility that x is a neutrino and concluded that a SU (2)l doublet neutrino was excluded 
by existing data. This was further reinforced by the subsequent improvement || in the experimental upper limit on 
the branching ratio, 
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BR(tt + -► fi+ + x) < 1.2 x 10~ 8 (95% C.L.), (2) 

versus the minimum value of ~ 2 x 10 -8 required in the doublet neutrino interpretation. However a sterile neutrino 
interpretation was found to be consistent, within strict limits on the mixing parameters (see also Ref. J?j), although 
this may still be in conflict with astrophysical and cosmological constraints ||. 

In Ref. H a solution was proposed based on the anomalous muon decay /i + — > e + + x, where x is taken to be a 
scalar boson of mass 103.9 MeV. However this implies too large a value for the energy released in the x decay and 
the required branching ratio is also constrained by the recent bound BR(/i+ — > e+ + x) < 5.7 x 10~ 4 (90% C.L.) [||. 
Thus it is necessary to add to the model 2 other scalar bosons into which x can cascade decay in order to dilute the 
energy This model appears viable but is somewhat baroque. 

In Rcf. Q, a supersymmetric solution was considered. The x particle was interpreted as a photino (or zino) and 
the anomalous pion decay 

7T+ -> fl+ + 7 (3) 



x The required mass is within 0.02% of the pion-muon mass difference. 
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was assumed to proceed via the R-parity violating operator L2Q1DI. (] The same operator then enables the photino 
to decay radiatively as 

7 7 + ^ (4) 

via a one-loop diagram with a d quark and d squark in the loop. However the expected peak at 17 MeV has not been 
reported in the new data Q on the energy spectrum of the anomalous events, so a 2-body decay for the x particle 
seems disfavoured. Therefore this model || may not be viable in its present form. We present below the necessary 
extension to produce a 3-body decay for such a light neutralino. n 



II. THE MODEL 



We consider the lightest neutralino in supersymmetry, x{> to be the hypothetical x particle, with mass m^o = 
33.9 MeV. This will also be the lightest supersymmetric particle (LSP) in our model. Since x is effectively stable on 
collider time-scales (this is quantified below) our model will experimentally look very similar to the MSSM. Now in 
a GUT-inspired MSSM, Mi — (5/3) tan 2 9^M2, and assuming this relation requires m^o > 32.3 GeV from current 
LEP data |L3]]. Thus in order to obtain a very light neutralino we must consider Mi and M% to be independent 
parameters. A small M2 implies at least one light chargino which is excluded by experiment, while a small Mi implies 
that the LSP will be dominantly bino. We will quantify this below and determine regions in the (Mi, M2, H, tan (3) 
parameter space consistent with all experimental limits. The solutions indeed turn out to be dominantly bino with a 
small higgsino contribution. 

Furthermore we invoke 2 non-zero R-parity violating operators. The pion decay 



(5) 

proceeds through the operator X' 211 L2QiD1 and the leading order Feynman diagrams are shown in Fig. |l|. The 
neutralino is assumed to decay as 



X? 



(6) 



through either \121L e L ^E^ or XisiL e L T E^. Note that this is the only kinematically accessible tree-level 3-body visible 
decay for such a light LSP. 




FIG. 1. Tree-level Feynman diagrams for pion decay via the operator LiQ\D\. 



In Fig. U we show the values of the branching ratio for 7r + — > /i + + x\ and lifetimes t-^o which are compatible 
with the KARMEN data In order to determine the required range of the couplings A' 211 , Ai{2,3}i in our model 

which are consistent with the solutions in Fig. || we must first determine the pion branching ratio in terms of the 
supersymmetric parameters. The partial width as computed from the diagrams in Fig. |l| is 



r(7r - mxi) 



87r(m d + m„) 2 
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(7) 



2 For reviews on R-parity violation see Ref. 



Such a decay was also proposed in Ref. llOl which invoked possible mixing between neutrinos and gauginos/higgsinos as 
the reason for neutralino instability rather than R-parity violating vertices. However to explain the KARMEN anomaly then 
requires the Higgs mixing term /1 H 1 H2 in the superpotential to be unnaturally small, /i < 30 MeV. Moreover this scenario 
implies a MeV mass v T which is definitively ruled out by cosmological and astrophysical arguments ]ll| , p^ |. 
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FIG. 2. The solutions to the KARMEN anomaly in terms of the anomalous pion branching ratio and the a;-particle lifetime. 
The hashed area denotes the experimental upper bound (0) . 



where m» and wip denote the pion and the muon masses, M~ - j, denote the corresponding scalar fermion masses, 
m u , rrid are the first generation current quark masses, and /„. is the charged pion decay constant. The constants 
A-e,u,d refer to the neutralino coupling and are given in Table | both for the general case and for the limiting cases of 
either a pure bino or photino neutralino. The phase space factor is given by p cm — ([m% — (r 

m f) 2 ]) 1 ' 2 /P ra i)- In the Appendix we give some details of how this result is obtained. 



The branching ratio for the anomalous pion decay, assuming the decay 7r + — > to be dominant, is given by 



BR(tt -> MX?) 
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(8) 



(9) 



M f , 



To obtain a numerical estimate, we have assumed in Eq. (|9j) that the scalar fermions are mass degenerate, M~ - j 
and that the neutralino is pure bino. In the last line we have quoted the experimental bound (||), shown as a hashed 
area m Fig. |. This bound can be satisfied by a small coupling and/or a large sfermion mass. It can also be satisfied 
by a fine-tuned cancellation between different diagrams for distinct sfermion masses, but we disregard this possibility. 
The last inequality (|^) can be translated into an upper bound on \' 2 n'- 



4 This corrects the result given earlier 



TABLE I. Neutralino coupling coefficients for the pion decay. 
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TABLE II. Coefficients for the neutralino decay 



Coefficient 




General formula 




Pure photino Pure bino 


B 1 
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< 6-8 x 10- ( J ^y . (10) 
In the limit where 3> ti^o, the neutralino decay rate for the operator XijiLiLjEf, j = 2, 3, is given by p[ 
rf .-.o :ir+f .^ A ^ m l? p! , , Sift B lB , B 2 B; 
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(11) 



where Mg L ,g R pih denote the scalar lepton masses, and -Bi.2,3 are the relevant Xiff couplings given in Table ||. In the 
second equation we have again assumed a pure bino LSP and degenerate scalar fermions. As a numerical estimate 



for the lifetime of the bino LSP with m^o = 33.9 MeV we obtain 

2 



MmJ \150GeV 
< 4.78 x 10 2 s, (13) 

where the last inequality is obtained by using the bound (||) and the set of solutions shown in Fig. ||. The resulting 
bound on the coupling is 

A 1{2 , 3}1 >1.66xl0-3(_!|_) 2 . (14 ) 

Given the perturbative upper bound, \jk < V^Lr, and the lower bound on the sfermion mass from LEP 2, M; > 
100 GeV, we also have a lower limit on the lifetime of Tbi no > 2.6 x 10~ 4 s. Thus there is a solution range of 6 orders 
of magnitude in lifetime or 3 orders of magnitude in coupling. For these lifetimes the LSP is stable on collider physics 
time scales. 

We now have all the ingredients to fix the model parameters. In our model, each point alongthe curves in Fig. |^ 
corresponds to a specific anomalous pion branching ratio (^) and a specific neutralino lifetime ( |l2| ) . If we assume the 
scalar fermions are mass degenerate, we can translate this into specific values of A 211 and A^^i for a fixed sfermion 
mass. This set of solutions in the R-parity violating parameter space is shown in Fig. [5] for Mj = 150 GeV (solid line), 

Mi = 300 GeV (dashed line), and M; = 1000 GeV (dot-dashed line). The hashed lines at A, A' = y/In denote the 
perturbative limit. For large scalar fermion masses (> 1 TeV) we quickly run out of room for perturbative solutions. 
The solutions above and to the left of the stars are excluded by the inequalities (|l^, |l4|) . 



III. CONSTRAINTS ON THE R-PARITY VIOLATING COUPLINGS 

The R-parity violating couplings we have introduced violate lepton number and are thus constrained by laboratory 
experiments. The best bounds at the 2a level have been been summarized as |l5| 

X' 2U < 0.059 ' dR 



\100 GeV J ' 

Am < 0.049 ( 10 Q f g eV ) => Tbino > 0.24 s, (15) 

A131 < 0.062 f .JfiO r bino > 0.15 s. 
V 100 GeV / 
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FIG. 3. Solutions to the KARMEN anomaly in terms of the R-parity violating couplings \2nL2Q1D1 and Ai{2,3}i> for 
different (assumed degenerate) sfermion masses. The hashed lines indicate upper limits on the couplings from perturbativity. 
The stars and diamonds (squares) give the upper limits on the couplings A 2 n and A121 (A131), respectively. Solutions above 
and to the left of the stars are excluded, as are solutions below and to the right of the squares (diamonds). 



The bound on A 211 is from measurements of R^ = T(ir — ► eiy)/T(ir — > fiu) fig] , the bound on Am is from charged- 
current universality Jl6|], while the bound on A131 is from a measurement of R T = T(t — > evv)/T(r — > fivv) fujl . The 
above bound on A 211 is weaker than the bound (fl0|), and we do not consider it further. In Fig. g the above bounds 
on the coupling Am and Am forbid solutions to the right of the diamonds and squares, respectively. We are thus 
left with a range of solutions of about 2 orders of magnitude in A' 211 and Au2,3}i- This corresponds to 4 orders of 
magnitude in the pion branching ratio and in the LSP lifetime, respectively. In the last 2 equations we have translated 
the upper bound on A!{ 2 .3}i into a lower bound on the lifetime using Eq. @, to be compared with the upper bound 
(fl3|). Note that these bounds are independent of the sfermion mass. 

Besides bounds on individual couplings, we must also consider bounds on the product of the couplings A^Am or 
A211A131. In the first case, we can get an additional contribution to pion decay ir + — > jl + — > e + v e which changes the 
prediction for R n 



r> _ R SM 
tin — 



TO 7r^'21lAl21 



2y/2G F M? L m e {m u + m d ) 



(16) 



As the corresponding Feynman diagram has a different structure from the t-channcl squark exchange which gives 
the bound on A'211 in Eq. ( |l5| ) we get a much stricter bound on the product of the couplings than on either of the 
couplings individually. This leads to the following bound at the 2a level 



This means that in the case of Am (as opposed to Am) the maximum scalar fermion mass which will solve the 
KARMEN anomaly in our model is 450 GeV. 

The couplings A 211 and Am violate muon and tau lepton number, respectively, and can thus lead to the decay 
t — * fvy. The experimental bound has recently been improved [fj"7f , 

BR(r -> ^7) < 1.0 x 10~ 6 (90% C.L.), (18) 

but is still 4 orders of magnitude weaker than the experimental upper bound on BR(/i — > ej). Therefore a bound 
on a product of couplings which yield the decay r — > /17 via a one-loop penguin diagram, e.g. Am Am, must be 2 
orders of magnitude weaker than the corresponding bound on the couplings which give fi — > ej, i.e. one would expect 
Ai 2 iAi3i < 0(1O~ 2 ) In our model, the couplings A 2U Am only contribute to the decay r — > /17 at the 2-loop 

level and the bound is thus significantly weaker than 0(1O~ 2 ). (The decay t — * fj.ee is similarly suppressed.) This is 
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significantly weaker than the bound (15) so we have no new bounds on the product X^n^iai- Furthermore since the 
bound (|l^) on A' 211 is so restrictive in our model, we need not worry about the model dependent bounds from flavour 
changing neutral currents |l9| . 

In Refs. po|-p2| severe cosmological bounds were derived on all R-parity violating couplings from considerations of 
GUT-scale lepto/baryogenesis in the early universe: 



A, A', A" < 5 x 1CT 7 i^j^) ■ (19) 

Subsequently it was shown that it is sufficient for just one lepton-flavour to satisfy this bound ]2l]^^ |. In Fig. |^ we 
can see that for our model both couplings violate the bound (|l9|). For the case (A' 211 , Am) we must therefore demand 
that either all electron number violating couplings or all tau number violating couplings satisfy Eq. (|l9|), while for 
the case (A 211 , A131), we must demand that all electron number violating couplings satisfy Eq. (JlQj) . Alternatively, 
baryogenesis could plausibly occur at the electroweak scale, in which case the bounds (|l9|) do not apply. 



IV. EXPERIMENTAL CONSTRAINTS ON A LIGHT NEUTRALINO 



We now summarize relevant experimental constraints on a light neutralino LSP and show that these are satisfied 
in regions of (Mi, M 2 , fi, tan/3) parameter space for a dominantly bino Xi with a small higgsino contribution. In our 
model, M\ and M 2 are not related by the supersymmetric grand unified relation and we treat them as separate free 
parameters. 



A. Bounds from e + e — > vv^ 

The process e + e _ — > 1/F7 can be measured in electron-positron collisions by detecting the photon and the missing 
energy due to the neutrinos |2^]. As the lightest neutralino in our model is long-lived on the time scale of collider 
experiments, the process e + e~ — > XiXil wlu gi ye the same experimental signature. 

The cross section for this latter process has been calculated [Q for the case of a pure photino and can be easily 
extended to the pure bino case we are considering here. The cross section is shown as a function of the centre-of-mass 
energy in Fig. ^ and is rather low. The expected number of events for a number of different experiments are given in 
Table III assuming a scalar fermion mass of Mj = 150GeV. We have used the same cuts on the energy and angle of 



the photon as in Ref. |24| ]. 

As can be seen from Table [II, no limits on this process can be set by LEP, as the expected number of events is 
much less than one. The recent results from OPAL give 138 observed events (with a statistical error of ±11.9) 
against the standard model expectation of 141.1 ± 1.1 events from e + e~ — > vDj. There is also an expected non-physics 
background of 2.3 ± 1.1 events. Thus there is no evidence for any excess. 

With the higher luminosities expected at the B-factories KEK-B and BaBar, a few events may be expected. The 
Standard Model (SM) cross section at this energy is 2.3 fb, corresponding to 230 ± 15 events at KEK-B and 70 ± 8 



TABLE III. Cross-sections for the production of XiX?7 a t e+e colliders for (an assumed common) sfermion mass 
Mr — 150 GeV and general expectations of the integrated luminosity. 



Experiment 


Integrated luminosity (pb 1 ) 


Energy 


Cross-section (fb) 


Number of events 
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6.65 
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0.04 
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0.04 
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FIG. 4. Cross-section for for the production of a purely bino neutralino with mass 33.9 MeV through e + e~ 



xlxh- 



events at BaBar. ^| The statistical uncertainty still exceeds the signal rate so we do not expect any sensitivity to a 
light neutralino. 

At the NLC we expect a substantially higher number of events. The SM cross section for the same cuts is about 
0.35 pb for 3 neutrinos p7| corresponding to 1.1 x 10 s events, with a small statistical error of 330 events. Thus this 
can provide a test of our model. 



B. Bounds from the invisible decay of the Z° 

In our model, m^o <C M z o/2, therefore the decay Z° — * XiXi is kinematically accessible, and the Xi can be 
considered to be effectively massless, just like a neutrino. The LSP decays outside the LEP detectors, thus the 
process Z° — > XiXi wn l contribute to the invisible width of the Z°. The current measurement of the invisible Z° 
width translated into the number of light neutrino species is ||2q| 



3.00 ±0.08, 



(20) 



so we must require that T(Z° — > XiXi) < 0-08 T(Z° — * uv), where the rhs refers to one neutrino species only. 

A pure bino LSP does not couple to the Z° at tree-level. The dominant contribution will thus come from the 
Higgsino admixtures of the LSP, Ni%, N±±, in the notation of Ref. pq . This enters with the fourth power in the 
decay rate Z° — > XiXi- The Higgsino has equal strength coupling to the Z° compared to a neutrino, thus yielding 
the constraint 



^/\N^T\N^ < 0.5 w (0.08) 1 / 4 . (21) 
We shall see below that it is straightforward to find regions which satisfy this in (M 1; M 2 , fi, tan/3) parameter space. 



5 This corresponds to one year of running based on the luminosities given in Ref. [E6|. 
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FIG. 5. Solutions in (Mi, M2, /J,, tan/3) parameter space giving a m x o = 33.9 MeV neutralino for /i = 300 GeV and 2 
representative values of tan/3. The width of the lines is 0.01 MeV. Below the hashed lines the chargino mass is less than 
150 GeV. The dotted lines have ApsuSY < 10~ 4 and the solid lines have ApsuSY < 5 x 10~ 4 . 



C. Solutions in MSSM parameter space 

It is important to establish whether it is indeed possible to have a neutralino LSP with m x o = 33.9 MeV within 
the MSSM. To this end we have scanned the MSSM parameter space (Mi, M2, fi, tan/3) with independent Mi, M2, 
for a neutralino in the mass range 



33.89 MeV < m^> < 33.91 MeV. 



(22) 



This leads to the neutralino iso-mass curves shown in Fig. [|, taking fi — 300 GeV and 2 representative values of tan /?. 
We have not been able to find any solutions with /i < 0. In order to obtain such a light neutralino some fine-tuning 
is required in the MSSM parameters, of about a few parts in 10 3 for tan j3 = 1 and a few parts in 10 2 for tan j3 = 8 
p9| . The fine-tuning is reduced for larger M2 and fi and small Mi because a light neutralino can then be generated 
by the see-saw mechanism; it is reduced for large values of tan/3 because in the limit (3 = tt/2 there is a zero mass 
eigenvalue for Mi « p9[ |. 

We have checked that the Higgsino contribution always satisfies the bound (|l]). In order to avoid an observable 
light chargino we require m ± > 150 GeV, which eliminates the region below the hashed lines in Fig. rafor the specified 
values of tan j3. 

The LSP is indeed dominantly bino along the solution curves in Fig. ||. The second lightest neutralino, £2) i s 
dominantly wino for M2 < 300 GeV, while for larger values it is mainly higgsino. For M2 ^ 110 GeV, m^o > 100 GeV, 
and for M 2 > 235 GeV, m^o > 200 GeV. 



D. Bounds from oblique electroweak radiative corrections 



Any new physics which couples to the SM can give contributions to the electroweak precision observables via 
radiative corrections. The effect of the new physics on vacuum polarization diagrams, the so called oblique radiative 
corrections, is usually parameterised using either the S, T, and U parameters of Ref. |30[| or the ei, £2 and €3 parameters 
of Ref. plj . The calculation of these parameters is based on an expansion in q 2 /M^ ev! , where q 2 is the momentum 
scale of the gauge boson propagator, typically M§ or smaller and M 2 ew is the scale of the new physics, assumed to be 
well above M\ . If however there are new light particles in the spectrum, as in the present case, these approximations 
are typically insufficient and one must in general also calculate the box or vertex corrections p2| . This full calculation 
is however beyond the scope of our analysis and is not attempted here. 
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There is one exception however and that is the ratio of the charged to neutral current neutrino-electron/muon 
scattering events — the p-parameter. This is defined at q 2 = and the expansion is thus trivial. In the following we 
limit ourselves to a calculation of the contribution to the p parameter from the full set of charginos and neutralinos. 
The radiative correction to the p parameter, Ap, is given by the W and Z self energies with zero momentum flow [] 

Ml 1 > 

The dominant SM contributions to Ap arise from the (heavy) top quark and the Higgs boson. Assuming the mass 
of the latter to be Mjj = Mz, and subtracting the SM contributions we are left with the experimental limit on new 
physics at the 2er level of |26| 

- 3.7 x 1(T 3 < Ap now < 1.1 x 1(T 3 . (24) 

We have calculated the contribution to Ap, which we denote ApgusY, from all charginos and neutralinos for given 
parameter points (Mi, M2, p,tan/?). We find full agreement with the results given in Ref. p3| . We then determine 
ApsusY along the solution curves given in Fig. |^. The dotted lines indicate solutions for which ApsuSY < 1CT 4 , while 
the solid lines indicate solutions for which ApgusY < 5 x 10~ 4 . Thus there is no conflict at least with the experimental 
constraint on the p parameter. 

E. Cosmological and astrophysical constraints 

Massive particles are expected to come into thermal equilibrium in the early universe and their relic abundance is 
essentially the equilibrium value at 'freeze-out' when their self-annihilation rate drops below the Hubble expansion 
rate. For the neutralino under consideration, the self-annihilation cross-section is (s-wave) suppressed p4| so the 
surviving abundance is rather high: 

(^) « 1.2 x 10- GeV (33^)" (t^)' - 0*) 

This energy density will be subsequently released when the neutralinos decay and this has the potential to disrupt 
standard cosmology, in particular primordial nucleosynthesis p"l"[ . Specifically, since the neutralinos will be non- 
relativistic during nucleosynthesis, the Hubble expansion rate will be speeded up, while the electromagnetic energy 
generated through the subsequent decays will dilute the nucleon-to-photon ratio, resulting in an increased helium-4 
abundance [E35). The decay electrons will also Compton scatter the thermal background photons to energies high 



enough to directly alter the abundance of e.g. deuterium through photodissociation |36|. The observationally inferred 
primordial abundances thus enable stringent bounds to be placed on the relic abundance of the decaying particle 
as a function of its lifetime. For the above abundance (|25|), the decay lifetime is required to be less than a few 
thousand seconds in order that the primordial D/H ratio is not reduced below its conservative lower limit of 10 -5 , 
and further required to be less than a few hundred seconds in order that the primordial 4 He mass fraction not exceed 
its conservative upper limit of 25% p5|] . Thus the cosmological lifetime bound is essentially the same as the one 
derived earlier ( |13| ) from experimental considerations. 

Very light neutralinos can also be produced through nucleon bremsstrahlung and e + e~ annihilation in supernovae 
such as SN 1987a. Since the squark/selectron masses are now restricted to be above raw |@, the neutralinos cannot 
be trapped in the supernova core by scatterings on nuclei or electrons, so will escape freely. The energy lost through 
this process can be comparable to the neutrino luminosity so may result in significant shortening of the i7 e burst. The 
neutralino luminosity can be decreased by increasing the sfermion mass but it has been shown that consistency with 
observations of SN 1987A is not possible for any sfermion mass less than 0(1) TeV jj7|. This constraint is evaded if 
the neutralino is unstable due to R-parity violation, as in the present case. However there are then further constraints 
on the energy released in the decays. Given the experimental upper bound (Ir3) as well as the cosmological upper 



bound on the lifetime, the decays would have occurred within the progenitor star. Moreover the lower bound ( |15[ ) on 
the lifetime implies that the neutralinos cannot have decayed within the supernova core. The electromagnetic energy 
released in the decays would have been thermalised leading to distortions of the lightcurve. However the neutralinos 
under consideration here have a mass which is of the same order as the core temperature [I4I so one must reconsider 
their production rate in order to quantify this potentially important constraint. 



3 See for example the third paper in Ref. p(l for a derivation of this result. 
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V. OTHER IMPLICATIONS FOR R-PARITY VIOLATING PHENOMENOLOGY 



A. The HERA high-Q 2 anomaly 

In 1997, the HERA collaborations reported an anomaly in their high Q 2 data @j whose most likely explanation was 
in terms of R-parity violation p9| |. However, this required a significant L e QiD^ operator. While this is not excluded 
by our model it is a completely distinct possibility. Together the two operators might contribute to the decay fj, — > ej. 
However this does not lead to a significant new bound since the bound on the coupling X' 211 is already so strict Jl8| . 



B. Neutrino masses 



The trilinear R-parity violating couplings we have introduced also generate Majorana masses for neutrinos through 
one-loop self-energy diagrams p9: 



^—^ 167T 2 m?~ 

k,a,b f 



(™t R )abm fk , (26) 



where N c is a colour factor, and (™|,fl)afc is the left-right mixing term in the sfermion sector. The question then 
arises whether our model can account also for the SuperKamiokande data suggesting oscillation of atmospheric 
into v T jil). The results indicate that the neutrinos mix almost maximally and that they are nearly mass-degenerate, 
5 m 2 ~ 10 -3 — 10 -2 eV 2 . If neutrino masses are hierarchical then the natural interpretation is that one of the neutrinos 
(presumably the v T ) has a mass of O(0. 1) eV (although the possibility of a close mass-degeneracy for a heavier v^—v^ 
pair is not excluded). Now the sfermion left-right mixing is not well determined, however within a given framework, 
e.g. supergravity-inspired models, approximate relations such as (rn 2 LR ) aa ~ m /a m / arise. Thus A 211 , the coupling 
responsible in our model for the pion decaying to the neutralino, generates a mass 

m MM « 1.5 x 10- 7 eV ( ' ( -^—) , (27) 



,io- 4 y V 150GeV , 

which is too small to be of phenomenological interest. The couplings A!{ 2 ,3}i responsible for neutralino decay also 
generate rather small masses: 

»w- ^.v^HisSb)". <»> 

Thus the absolute scale of the masses seems too low to explain the data on atmospheric neutrinos. However if other 
R-parity violating couplings are also present, it may well be possible to generate a neutrino mass pattern consistent 
with the observations, both of atmospheric and solar neutrinos Q . 



VI. FUTURE TESTS 



Experimentally, our model largely looks like the MSSM with non-universal gaugino masses and with a very light 
LSP. Thus most future tests of the MSSM also apply to our model, e.g. chargino pair production. A specific test 
would be to identify a very light LSP for example via neutralino pair production p 7 ^ . At an e + e _ collider one can 
study the process 

e + + e-^x%+xl (29) 

where x 2 subsequently decays visibly Q]. In Fig. ^ we show the cross section evaluated along our MSSM solution 
curves for both LEP2 (-y/s = 200 GeV) and the NLC (-y/s = 500 GeV). This should be directly observable, provided 
it is kinematically accessible, i.e. m-^o < y/s. 

Experimentally the main difference between R-parity violation with a long-lived neutralino LSP and the MSSM is 
the possibility of resonant sparticle production. The value of A^n ^ 10~ 4 is too small for the observation of resonant 
slcpton production at hadron colliders However values of A!{ 2 ,3}i > 10~ 3 should allow a test for resonant second 
or third generation sneutrino production at e + e~-colliders for masses upto ~ y/s p6| . One can also test for the first 
generation via the mechanism described in Ref. p7|. 
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FIG. 6. Cross-Sections for e + e~ — > X2X1 f° r the solutions in Fig. ^ The solid lines correspond to 0.1 pb< a <1 pb, the 
dashed lines to 10 fb< a <0.1 pb, the dotted lines to 1 fb< a <10 fb, and the dot-dashed lines to a <1 fb. 



A further upgrade of the KARMEN detector may allow a better resolution of the decay of the x particle, in 
particular the angular distribution of the decay products. For reference we show in Fig. [7] the differential decay rate 
of the LSP in our model as a function of the angle between the two final state electrons. 



VII. CONCLUSIONS 



The KARMEN time anomaly is particularly intriguing because contrary to several other reported 3 — 4cr effects 
in the literature, its significance has not diminished with improved statistics, nor has it been explained away as a 
systematic effect. In fact the anomaly persists in the KARMEN-2 data, which has a much reduced background [||, 
with the same characteristics as in the KARMEN- 1 data 0. It would appear that there is no independent experiment 
which has the sensitivity to reproduce this result. In particular although the LSND experiment also studies pions and 
muons decaying at rest, it lacks the distinctive time-structure of the beam in the KARMEN experiment necessary to 
isolate the anomaly. Since KARMEN-2 acquires only of O(10) anomaly events per year of running, it is clear that a 
definitive resolution of the problem will have to await an upgraded detector with tracking capability. 

Phenomcnological models for the anomaly as due to the production and decay of a new particle are very tightly 
constrained. The only viable proposals at present are a singlet neutrino decaying through its large mixing with the 
v T , or a neutralino decaying through violation of R-parity Q] which we have extended and investigated in detail. 
An important lesson from our investigation is that contrary to popular belief a neutralino lighter than even the pion 
is not excluded by present accelerator data unless a GUT relation between gaugino masses is assumed. Whether the 
KARMEN anomaly is indeed the first evidence for such a particle is a matter for future experiments to decide. 



Note Added: While this manuscript was under review, the E815 (NuTeV) experiment at Fermilab reported a search 
for a 33.9 MeV neutral partic le produced in p ion decay decaying to a partially electomagnetic state such as e + e~v 
or -fv (J. A. Formaggio et at, hep-ex/9912062 ). No evidence was found for such a particle but the lifetimes probed 
(~ 1CU 9 — 1CU 3 s) are much smaller than the lower limits ( |l5| ) on the neutralino lifetime in our model so there are no 
implications. We note however that the exclusion of such short lifetimes is relevant in the context of the constraints 
from SN 1987a on the decaying particle hypothesis ||[||- These constraints have been investigat ed further in two 
other recent reports (I. Goldman, R. Mohapatra and S. Nussinov, hcp-ph /99 1 2464 M. Kachelriess, ^icp-ph/0001160 ) 
which conclude that our model is excluded by the observations of SN 1987a. We reserve judgement on this issue for 
the reasons mentioned earlier. 
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FIG. 7. Decay rate of the LSP versus the angle between the final state electrons. 
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APPENDIX A: CALCULATION OF THE PION DECAY RATE 

The rate of the pion decay it — > /i + Xi can be calculated using chiral perturbation theory. To do this we need an 
effective Lagrangian for the 4-fermion interaction u, d, Xi and /z + with the sfermion degrees of freedom integrated 
out (analogous to using the Fermi Lagrangian in the SM calculation). This gives 

c = >Wni ( ™ o) W) ' (A1) 

\ M u d / 

where we have also Fierz-reordered the Lagrangian (and neglected some tensor-tensor interaction terms which cannot 
contribute to the pion decay rate) . The matrix element of the axial- vector current between the pion and the vacuum 
is 

(0\r 5a (x)\TT b (p)) = -ip^Ud^e-^, (A2) 
where a and b are isospin indices. Using Eq. ( |A2| ) we obtain 

(0\u^ l5 d\n-) = -iVV/^e-*". (A3) 
Contracting this with the pion 4-momentum and using the Dirac equation for the up and down quarks yields 

(0^ 75 dK-) = + i ^< e ~[ PX , (A4) 
(m d + m u ) 

where md, m u are thus the current quark masses. The amplitude for the decay (|^) is then, 

^2\\^ rr ^n I A e A u Ad 

{m d + m u ) Vm? ~ 2M| " 2M 



A = J, I T72 - - T^h ) (uMPnvM) (2ir) 4 S (p - Pl - p 2 ) ■ (A5) 



This gives the decay rate quoted in Eq. (^). The additional contribution to the decay rate 7r H 
Section IB can be calculated in the same way. 
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